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AERODYNAMIC DAMPING AT MACH NUMBERS OF 1.3 AND 1.6 OF 
A CONTROL SURFACE ON A TWO -DIMENSIONAL WING 
BY THE FREE -OSC ILLATION METHOD 
By W. J. Tuovila arid Robert W. Hess 


SUMMARY 


Tests have been made at two supersonic speeds to obtain experimen- 
tally the aerodynamic damping characteristics of a control surface on a 
two-dimensional wing. The control surface had a chord of I .67 inches 
(l/3 of the wing chord) and a span of 7.25 inches and was supplied in 
three materials (steel, aluminum, and magnesium) having different mass, 
inertia, and stiffness properties. Two wing sections were tested, one 
being of 65A004 section and the other a 5 -percent-thick hexagonal sec- 
tion. The test results are compared with results calculated by two- and 
three-dimensional oscillating air-force theories. At a Mach number of 
1 . 6 , both theories are in fairly good agreement with the experimental 
results. At a Mach number of 1.3 , both theories predict negative 
(unstable) damping, whereas the tests indicate that the damp in g is 
slightly positive (stable). The in-phase or aerodynamic stiffness coef- 
ficients predicted by both theories are slightly higher than the experi- 
• mentally determined coefficients. 


INTRODUCTION 


Theoretical studies have indicated that at low supersonic speeds 
control surfaces with a single degree of torsional freedom can encounter 
unstable aerodynamic damping at some values of reduced frequencies.- 
Since existing theories do not account for many flow effects which may 
influence the problem, tests were made to obtain some exper ime nt ally 
determined aerodynamic damping coefficients for comparison wittiitoieo- 
retical values. Aerodynamic in-phase or stiffness coefficients and out- 
of -phase or damping coefficients were determined for a l/ 3 -chord control 
surface attached to a two-dimensional wing at zero angle of attack. 

Wings with hexagonal and. 65A004 section shape were used. The tests were 
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made at Mach n umb ers of 1.3 RT1 d 1.6 over a reduced frequency range from 
0.029 to 0.074. This paper presents the test results and compares them 
■with results calculated using two- and three-dimensional theories for 
osc ill ating air forces. The test results are also compared with the 
results of some damping tests made on a control surface attached to a 
triangular wing (ref. l). 


SYMBOLS 


-A-a 


aspect ratio of control surface, l a /2b a 

^a 


semichord of control surface, ft 

8 


control-surface deflection, radians 

fo 


natural frequency of rotation of control surface about 
hinge line at zero airspeed, cps 

ft 


natural frequency of rotation of control surface about 
hinge line at test Mach number, cps 

So 


damping coefficient associated with f Q 

st 


damping coefficient associated with ft 

^a 


length of control surface, ft 

ma 


mass of control surface, slugs/ft of span 



in-phase aerodynamic coefficient per foot of span 

% 


out-of -phase or damping coefficient per foot of span 

M 


Mach number 

V 


airspeed, fps 

P 


• air density, slugs/ft^ 

.“to 


2rtf 0 



atft 


k 


reduced frequency b a p%JV 
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xO 


C b£ 

la 

r 2 
r a 


spring constant, ft-lb/radian 
-2kN6 

mass moment of inertia about binge line , slug-ft 2 /ft of span 



“a/^a 2 


Subscripts : 
a 
o 
t 


refers to control surface 

refers to conditions at zero wind velocity 

refers to conditions at test Mach number 


MODELS AND TEST METHODS 


Wing, control surface, and binge details are given in figure 1. 
Control surfaces made of steel, aluminum, and magnesium were tested on 
two steel wing models which differed only in section. One wing model 
had a 65A004- section and the other had a 5-percent-thick hexagonal sec- 
tion. Each wing had a 5-inch chord and spanned the tunnel test section 
with one end clamped in the sidewall and the other end pinned in the 
sidewall. The control-surface chord was l/j of the wing chord. Steel 
hinges of various stiffnesses were used to attach the control surfaces 
to the wings at three points. There was a gap of about 0.02 inch between 
the wing and the control surface. Table 1 lists some of the physical 
parameters of the models. The masses and inertias were determined exper- 
imentally and include the contribution of the hinges. 

The tests were made in the 9- by 18-inch Langley supersonic flutter 
apparatus which is an intermittent-flow blow-down t unn el operated at 
atmospheric stagnation pressure. . The testing technique used was first 
to obtain "no-wind" damping decrements with the wing in the test ing con- 
figuration by flicking the control surface. The control surface was 
then deflected, the tunnel was brought up to speed, and the control sur- 
face was released and a "wind-on" damping decr emen t was obt ain ed - The 
air flow was then stopped' and the process was repeated using hinges of 
different stiffness. 
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Tire initial amplitude of both, the "no-wind" and the "wind-on" 
oscillations was not controlled precisely. It was judged by eye to 

range from about ±1° to ±2^°, the larger amplitudes occurring at the 

lowest frequencies. 

The system for deflecting the aileron is illustrated in figure 1 
Q-nfl consisted of a wire with an eye on the end which was inserted through 
a Rmsii hole at the trailing edge of the control surface. A straight 
release wire was then inserted through the eye of the cocking wire. The 
control surface was cocked by pulling the coc kin g wire until the desired 
deflection was obtained. The control surface was released by pulling the 
release wire out of the eye of the cocking wire. 

Damping decrements were obtained from a strain gage glued to a thin 
metal strip fastened to the wing and control surface. This metal strip 
followed the control -surface motion and the strain-gage output was 
amplified and fed into a recording oscillograph. 


SEDUCTION OF DATA 


The experimental decay decrements^ were reduced to average total 
supersonic aerodynamic coefficients N^ and Ng as was done in refer- 
ence 2 for subsonic flow. All damping terms are assumed proportional 
to amplitude and in phase with velocity. The following equation of 
equilibrium. 


I a 6 + k a (l + igq) 6 = -4pb a 2 V 2 k 2 8 ^ + iNgj 


leads to the following results for the in-phase component, 

21 


N 5 = pr a 


1- P 
w 


and, for the out-of -phase or damping component. 


% = nr a 


“ So 


(1) 


( 2 ) 


(3) 


The details of the analysis are given in the appendix. 
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It may be noted that the damping component is not obtained from 
just the difference in the damping coefficients of the "wind-on" and 
"no-vind" decrements. Instead, the "no-wind" damping coefficient is 
reduced by the factor which accounts for the difference in 

the structural damping coefficient due to the differ enc e in frequ en cy 
between "wind-off" and "wind -on" conditions. It is of interest to note 
that at M = 1.3 the "no-wind" damping coefficient t go was usually 


larger than the "wind-on" damping coefficient 

K/“t ) 2 

positive. 


gh but the factor 
made the aerodynam i c damping coefficient Eg slightly 


The experimentally determined E^ and Eg are compared with two- 

end three-dimensional air-force coefficients obtained from references 3 
and t . For comparison with the results obtained in reference 1, the 
damping coefficient Eg is expressed in stability notation using viscous 

type damping terms as follows : 


^ = - 2 “6 
0 V 

RESULTS AED DISCUSSION 

Presentation of Data and Comparison With Theory 


oo 


The control surfaces were attached to two-dimensional wings set at 
zero angle of attack. The aerodynamic in-phase and damping coefficients 
were obtained from the decay records and frequencies obtained in both 
still air and at the test Mach numbers of 1.3 and 1.6 and the data are 
presented in table 2. Sample "wind-off " and "wind-on" decrements are 
shown in figures 2(a) and 2(b). The hinge axis was so near the leading 
edge of the control surface that it wa s assumed to be there. The aero- 
dynamic damping coefficients Eg are presented in figure 3 and the in- 

phase coefficients N^ are presented in figure 4. The aerodynamic 

coefficients are plotted against the reduced frequency, based on the 
control-surface semichord. 

The experimental results are compared with the two-dimensional 
theory of reference 3 by assuming the control surface to be a wing 
oscillating about its leading edge and with the three-dimensional theory 
of reference k, assuming a sealed gap between the wing and the control 
surface. The theoretical results are also plotted on figures 3 and 4. 
Both theories predict negative aerodynamic damping at M = 1.3} however, 
the three-dimensional theory predicts only about l/2 the damping of the 



6 


HACA EM L56A26a 


too -dimensional -theory. The experimental aerodynamic damping at M = 1.3 
is slightly positive and both theories approach it as k increases. At 
M = 1.6 both theories are in good agreement with the experimental damping 
results, the three-dimensional theory giving slightly higher values than 
the two-dimensional theory. ■ 

The experimental in-phase aerodynamic coefficients presented 

in figure 4 are fairly consistent and both theories predict the trends 
well. The two-dimensional theory gives slightly higher values than the 
three-dimensional theory does and both theories yield values that are 
higher than the experimental. 

It appears that linearized flow theory, when applied to flow around 
trailing-edge control surfaces, begins to break down at low Mach numbers 
in the neighborhood of 1.3 or less. Adding an aspect-ratio correction 
to the two-dimensional -flow theory improves the results; however, some 
basic differences between the actual and the idealized flow appears to 
affect the results . Wing thickness, boundary layer, and the gap between 
the wing and control surface are some factors whose effects are not 
included in the theory. Also, the experimental results were obtained 
from decaying oscillations, whereas the theory assumes constant-amplitude 
oscillations. At M = 1.6 the theory seems to compensate for these 
effects and the agreement is good. 


Comparison With Control-Surface Data for a Triangular Wing 

The results of the present tests are compared in figure 5 with those 
of reference 1 through the Mach number range. Results for an amplitude 
of ±3° at a maximum k value of 0.03 from reference 4 are compared with 
the results of the present tests for amplitudes of about i2° at k values 
of 0.045 • The damping coefficients are expressed in stability notation 
as Qhg- The difference in the present results and those of reference 1 

may be the result of differences in flow caused by the wings. It may be 
noted that in reference 1 the control surface is attached to an aspect- 
ratio-2 triangular wing and not to a two-dimensional wing. In refer- 
ence 1 the damping varied from a small degree of instability at M = 1.3 
to neutral stability at M = 1.9, whereas the present tests indicate 
slight stability at M = 1.3 and considerable stability at M = 1.6. 

The too- and three-dimensional theory results are also presented in 
figure 5. 


GENERAL OBSERVATIONS 


At M = 1.3 there is considerable scatter in the results but the 
damping coefficients in all but one case are positive. This scatter is 
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due to the sensitivity of the equation for Eg to small changes in 

measured damping between the "wind-off" and "wind-on" conditions when 
the aerodynamic damping is low. Ho flutter was observed during these 
tests which indicates that the total damping was positive and shows 
that the aerodynamic damping could have been, at most, only slightly 
negative since the structural damping was small. At M = 1.6, where 
the aerodynamic damping is higher, the scatter is considerably reduced. 

Any effects due to wing-profile or control-surface material is lost 
within the scatter of the results. 

The structural damping g Q was principally in the range 0.006 to 
0.01 with a few extreme values of gQ = 0.004 on the low end and 
g Q = 0.034 on the high end. This spread in the structural, damping 
coefficient is believed to be due to variations in the hinge clamping 
force. Also, the structural damping coefficient generally decreased 
with decrease in amplitude and some unusually large changes are noted 
in table 2(a) for pr 0 ^ = 650 . The damping coefficients recorded in 

table 2 (a) were measured near the maximum amplitude of oscillation. 

The aerodynamic dancing may also be affected by amplitude; however, 
since the present tests were made without amplitude control, no such 
effect can be determined. No appreciable amplitude effect is indicated 
in reference 1 at Mach numbers from 1.3 to 1.9 while reference 5 shows 
considerable effect for amplitudes up to ±5° at Mach numbers near 1.0. 

Wing bending motion may also affect the results by introducing a 
translation degree of freedom to the control surface. Although the wing 
motion was not measured, it is believed to have been very slight since 
the wing was clamped at one end and pinned at the other. As the control- 
surface frequency approached the wing resonant frequency, the wing ampli- 
tude would increase rapidly and any bending effect should become evident. 
At M = 1.6 the HACA 65A004 wing with control surface pr & 2 = 378 

reached the wing resonant frequency at k = 0.069 and yielded essentially 
the same results as the hexagonal wing with control surface pr a ^ = 4-27 

where the control-surface frequency was 85 percent of the wing res onant 
frequency. The HACA 65A004 wing would have had, about 5 times the ampli- 
tude of the hexagonal wing at this k value which ind icates that the 
wing bending amplitude had no apparent effect on the damping results. 


CONCLUDING E1MARKS 


The results of the tests of a control surface attached to a two- 
dimensional wing at zero angle of attack ind icate that at a Mach number 
of 1.3, a slight amount of aerodynamic damping exists on the control 
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surface, whereas both two- and three-dimensional theories predict negative 
damping. At a Mach number of 1.6 the control surface has considerable 
aerodynamic damping which both two- and three-dimensional theories predict 
quite well. Both theories predict the trends of the in-phase aerodynamic 
coefficients, but they yield results which are slightly higher than 
experimental values. These results were obtained at reduced frequencies 
from 0.029 to 0.07^. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., Janua r y 9 , 195 6. 
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APPENDIX A 

Derivation of Aerodynamic Coefficients N^ and Ng 

The supersonic aerodynamic coefficients N^ and Ng are derived 

from the following equations of equilibrium, where the damping is assumed 
proportional to the displacement and in phase with the velocity. 

"Wind-on" condition (aerodynamic and structural) 


I a 5 + kt(l + ig t )s = 0 

(Al) 

"No-wind" condition (structural only) 


I a 5 + ka (l + igo) 5 = 0 

(A2) 


where 


itg = Logarithmic decrement 

then, 

kj.8 - k a S = Aerodynamic spring force 

and, 

k^g^S - k a go5 = Aerodynamic damping force 

Equation (a 4) implies that the structural damping force is 
independent of frequency. 

By definition 


Aerodynamic spring force 

fob -*a)B 


(A3) 

m 


(A3)' 


% 


_ Aerodynamic dancing force 
i<-pb a 2 V 2 k 2 6 

= - ^agp) 5 

4pb a 2 V 2 k 2 5 


(A6) 


and 
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for small values of damping 


H l*a> % 2 = 1%/la 

reduced frequency (by definition), 

k = b a o^ V 

Substituting equations (AJ) and (a 8) into (A5) and (A 6) 


n 5 = 


^(“t 2 - “b 2 ) 


^P b a V ‘ 


2xr2 b a 2f % 2 


1 - P) 

W 


V 2 


- _ Ia(o>t 2 gt - Ob 2 go) 

6 - "2Z2 


72 


_ *a 


4pb. ^ 


a i— 


St “ So 


©' 


finally, substituting 


P-r a 


2 


_ma_ x 

. 2 

4pb a 


la 

2 

“a^a 



% = 


1 - f2a 
K 


% = ^a C 


- So t 


(AT) 

(A8) 

(A9) 

(A10) 

(All) 

(A12) 
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TABLE 1.- SOME CONTROL-SURFACE PHYSICAL PARAMETERS 


Wing section 

NACA 65A004 

Hexagonal 

Control-surface 

material 

Steel 

Aluminum 

Magnesium 

Magnesium 

•••••• 

0.606 

0.606 

0.606 

0.600 

b a •••••• 

O.O 696 

0.0692 

0.070 

0.0679 

m a •••••• 

0.0lh5 

0.00593 

0.00357 

O.OO679 

la •••••• 

5-5 X 10“5 

2.29 X 10“5 

1.24 X 10“5 

2.39 X 10“5 

r 2 

x a •••••• 

0.782 

0.806 

0.71 

O.766 

^M-1.6 * * * • 

1133 

469 

276 

559 


p = 0.00090 at M = 1.3 and V = 1^30. 
p = 0.00066 at M = 1.6 and V = 1760. 


The first natural wing frequency for the NACA 65AOC& wing was 
about 260 cps and for the 5 -per cent hexagonal wing it was about 
300 cps. 
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TABLE 2 .- CONTROL-SURFACE DATA 
(a) M = 1.3 



f t 

So 

®t 

k 

% 

% 

"CHS 



Hexagonal wing; 

^r a 2 = 

313 



68 

182 

0.0115 

0.0023 

0.0550 

269 

0.22 

0.024 

70 . 

182 

.0105 

.0022 

.0550 

267 

• 19 

.021 

68 

184 

.010 

.007 

• 0555 

270 

1.75 

. 194 - 

72 

183 

.0095 

.011 

.0552 

264 

2.98 

.330 

89 

188 

.0085 

.0082 

.0568 

242 

1-97 

.224 

101 

193 

.0095 

.0092 

.0585 

228 

2.07 

.242 

101 

191 

.011 

.007 

.0578 

228 

1.25 

.144 

133 

215 

.019 

.012 

.065 

193 

1.47 

.191 

133 

213 

.034 

.015 

.0645 

191 

.53 

.068 

145 

222 

.013 

.012 

.0672 

180 

2.04 

.274 

146 

217 

.014 

.010 

.0658 

171 

1.16 

.153 

178 

240 

.018 

.0115 

.0725 

l 4 l 

.50 

.073 

180 

242 

.014 

.0074 

.073 

l 4 o 

-.09 

-.013 

185 

245 

.0165 

.014 

.074 

135 

1.44 

.213 



65 A 004 wing; 

|ir a 2 = 650 



66 

135 

O.OO85 

0.0068 

0.0419 

4-95 

3.12 

0.262 

66 

136 

.0070 

.0030 

.0422 

4-97 

.85 

.072 

81 

142 

.0109 

.0085 

.0440 

439 

3.25 

.286 

81 

142 

.0097 

.0075 

.0440 

439 

2.80 

.246 

81 

142 

.0091 

.0052 

.0440 

439 

1.43 

.126 

89 

146 

.010 

.0093 

.0453 

408 

3.64 

.330 

89 

147 

.0075 

.0070 

• 0457 

412 

2.80 

.256 

89 

145 

.0090 

.0050 

.045 

405 

l.o 4 

.094 

120 

164 

.010 

.0055 

.0508 

302 

.06 

.006 

120 

165 

.010 

.0072 

.0513 

306 

1.24 

.127 

158 

175 

a. 032 

.021 

.0544 

246 

.65 

.071 



la.018 

.021 



6.50 

.706 

132 

172 

a. 018 

.027 

.0535 

267 

io. 4 o 

1.112 



b .oi8 

.021 



6.50 

.695 

147 

186 

.006 

.010 

.0578 

244 

4.03 

.465 

160 

192 

a. 020 

.028 

.0595 

198 

9.10 

1.082 



b.0l6 

.017 



3.90 

.465 

160 

190 

a. 024 

■ .030 

.0589 

189 

8.50 

1.000 



4.016 

.014 



li 95 N 

.230 


a High amplitude. 
^Lower amplitude. 
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TABLE 2.- COMKOL-SURFACE DATA - Continued 
(b) M = 1.6 


fo 

ft 

So 

St 

k 

% 

n 6 

-Obg 



Hexagonal vingj ura^ 

= 427 



70 

160 

o.oi 4 

0.021 

0.0415 

345 

7.80 

0.647 

72 

163 

.011 

.020 

.0424 

344 

7.65 

.648 

71 

160 

.0107 

.023 

.0415 

343 

8.90 

• 739 

90 

170 

.on 4 

.021 

.0442 

307 

7.60 

.671 

90 

171 

.0103 

.0205 

.0445 

309 

7-51 

.669 

92 

170 

.0104 

.024 

.0442 

302 

8.90 

.786 

92 

170 

.010 

.0275 

.0442 

302 

10.50 

• 929 

101 

172 

.0103 

.0256 

.0446 

280 

9-40 

.838 

102 

176 

.0118 

.0259 

.0456 

284 

9 . 4 o 

.856 

148 

210 

.0097 

.0245 

.0547 

215 

8 . 4 o 

.920 

148 

201 

.010 

.024 

.0520 

196 

7-95 

.82 6 

l 46 

202 

.0115 

.023 

.0524 

204 

7-25 

• 759 

147 

202 

.0078 

.0188 

.0524 

201 

6.27 

.656 

189 

233 

.0060 

•0155 

.0606 

147 

4.95 

.600 

187 

231 

.0065 

.018 

.0606 

148 

5.85 

.710 

187 

233 

.0060 

.023 

.0606 

132 

8.15 

• 989 

215 

250 

.0058 

.019 

.065 

117 

6.36 

.827 

213 

251 

.0054 

.016 

.0654 

120 

5.16 

.675 

233 

265 

.0078 

.0214 

.069 

97 

6.57 

.907 

232 

260 

.0088 

.020 

.0676 

88 . 

5.55 

• 750 




MCA EM L56A26a 


15 


TABLE 2.- COUTROL-SURFACE BATA - Concluded 
Ob) M = 1.6 


f o 


So 


k 

n 5 

W 6 

-0u . 

h b 



65A004 wing; 

ur a 2 = 886 



52 

111 

O.OO78 

0.016 

0.0296 

693 

12.7 

0.752 

52 

112 

.0080 

.0154 

.0298 

696 

12.1 

.723 

53 

no 

.0073 

.0172 

.0294 

680 

13.7 

.805 

53 

113 

.OO67 

.0147 

.0301 

691 

H.7 

.705 

66 

120 

.0063 

.015 

.0321 

619 

n.6 

.745' 

66 

119 

.0066 

.0164 

.0318 

6l4 

12.8 

.815 

65 

120 

.0085 

.0168 

.0321 

626 

12.7 

.815 

65 

120 

.0083 

.016 

.0321 

626 

12.0 

.770 

80 

130 

.0076 

.0154 

.0347 

551 

n.i 

.770 

80 

129 

.0079 

.0164 

.0344 

545 

n.9 

.819 

80 

128 

.0060 

.0146 

.0341 

540 

10.9 

.744 

80 

129 

.0063 

.0158 

.0344 

545 

n .9 

.819 

90 

133 

.0062 

.0165 

.0354 

480 

12.1 

.856 

90 

133 

.0063 

.0151 

.0354 

480 

10.8 

.765 

90 

133 

.0055 

.0162 

.0354 

480 

12.1 

.856 

88 

131 

.0069 

.015 

.035 

487 

10.6 

.742 

152 

176 

.0047 

.0144 

.047 

225 

9-7 

.910 

153 

176 

.0040 

.0146 

.047 

217 

10.3 

.968 



65A004 wingj 

\ST a 2 = 378 



83 

176 

0.009 

0.026 

0.0467 

294 

9-13 

0.852 

83 

176 

.009 

.0223 

.0467 

294 

7.72 

.711 

105 

186 

.010 

.0224 

.0493 

258 

7.26 

.717 
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NA.CA EM L56A26a 



Figure 1.- Sketch of wing and control surface used in tests. 
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Control - surface decrement 





>.j |/60 second 

AAAAAAAAAAAAAA/VVWVW 


(a) "Wind off" decrement. 
Figure 2.- Sample decrement. 
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MCA RM L56A26a 
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AAAAAAA/WXAAAAAAAAA^ 


(b) "Wind on" decrement. 


Figure 2.- Concluded. 
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(b) M = 1.6. 


ynaralc coefficient TT^ vith reduced 
= 1.3 and 1.6. 
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S Present test results for k = .045 
□ Ref. I for wing at zero angle 
of attack and ±3° control 
surface amplitude 
O 3-D Theory (Ref. 4) 

O 2-D Theory (Ref.3) 


Unstable 





hQ- 




1.4 


Mach Number 


Figure 5*- Variation of damping coefficient with Mach number and 

comparison with results of reference 1 at k values of about 0.0?. 
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